Lagrangian statistics of light particles in turbulence 
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We study the Lagrangian velocity and acceleration statistics of light particles (micro-bubbles in 
water) in homogeneous isotropic turbulence. Micro-bubbles with a diameter db = 340 /im and Stokes 
number from 0.02 to 0.09 are dispersed in a turbulent water tunnel operated at Taylor-Reynolds 
numbers (ReA) ranging from 160 to 265. We reconstruct the bubble trajectories by employing three- 
dimensional particle tracking velocimetry (PTV). It is found that the probability density functions 
, (PDFs) of the micro-bubble acceleration show a highly non-Gaussian behavior with flatness values 

in the range 23-30. The acceleration flatness values show an increasing trend with ReA, consistent 
i-j^ " with previous experiments [ll and numerics @|. These acceleration PDFs show a higher intermit- 

r!jr , tency compared to tracers [3] and heavy particles in wind tunnel experiments. In addition, 

the micro-bubble acceleration autocorrelation function decorrelates slower with increasing ReA. We 
\^ " also compare our results with experiments in von Karman flows and point-particle direct numerical 

simulations with periodic boundary conditions. 

^3 ! I. INTRODUCTION 

^ ■ Multi-phase flows where the carrier fluid transports particles under turbulent conditions are ubiquitous. A 
thorough understanding of the dynamics of particles (light, neutral, or heavy) in turbulent flows is therefore crucial. 
In most of these flows, the particles have a finite size and their density is different from that of the carrier fluid. 
£>-j" Thus, the particle's dynamic behavior is expected to be different compared to neutral fluid tracers. The two relevant 
. dimensionless parameters are the density ratio ft — 3pf/(pj + 2p p ), where pj and p p are the fluid and particle density, 
and the Stokes number, which is the ratio of the particle's response time t p to the Kolmogorov time scale r^, defined 
as St= T p /r, ; = a 2 /ij3vT ril where a is the particle radius and v the kinematic viscosity of the carrier fluid. 
£N) ■ The Lagrangian approach is naturally suited to study particles in turbulence and has recently attracted much 
attention [seeQ]. Pioneering Lagrangian particle tracking experiments in fully developed turbulence used silicon 
OO . strip detectors to measure three-dimensional trajectories of tracer particles (/3=1) with high spatial and temporal 
resolution in a von Karman flow at high Taylor- Reynolds numbers, ReA up to 970 P, @,I3- The particle acceleration 
PDFs were found to be highly intermittent with flatness values around 55, and could be fitted with either stretched 
exponentials or log-normal distributions. The high intermittency of the fluid particle acceleration PDFs was also 
observed in numerical simulations @, I8l4l0| . The normalized acceleration PDFs showed a weak dependence on the 
Rca, and this was more prominently seen in the flatness values, which have been found to increase with ReA in both 
experiments [l[ and numerics Q. 

More recent experimental investigations have focused on studying particles with different density than the carrier 
fluid pd| - [T3 |. Heavy particles ((3 = 0) in turbulence were studied using water droplets in wind tunnel experiments 
at ReA=250 By following the particle motion with a moving camera, their trajectories were obtained in 

two-dimensions. It was observed that the normalized acceleration PDF was less intermittent than for tracers, with 
narrower tails. Numerical simulations [15| have confirmed that the acceleration PDF of heavy particles is indeed 
slightly narrower than that of fluid tracers. 

The dynamics of light particles in turbulence (/3=3) have been investigated both numerically and experimentally. 
DNS in the point-particle limit |9j showed very high intermittency in the PDFs of the individual forces acting on 
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bubbles. A direct comparison of the statistics of light, neutral, and heavy particles was done by Volk et al. [12j, using 
data from both experiments and point-particle DNS. In their experiments an extended laser Doppler velocimetry 
technique (extended LDV) was used to measure the particle velocity in a von Karman flow at high ReA =850. The 
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experimental PDFs of the normalized acceleration for light, neutral, and heavy particles did not reveal a visible 
difference within the experimental accuracy. In contrast, their numerical results showed that the acceleration PDF of 
light particles was more intermittent than that of tracers, and heavy particles showed less intermittency than tracers. 
They also found both in numerics and experiments, that the acceleration autocorrelation function of light particles 
decorrelates faster than those of neutrally buoyant and heavy particles. 

In this work, we present an experimental study of the Lagrangian dynamics of light particles in turbulence. A 
majority of the previous Lagrangian particle tracking experiments focused mainly on tracer particles. Furthermore, 
grid- gen erated turbulence experiments in wind tunnels have dealt with either heavy or neutrally buoyant particles 
[H, lil I 111] . Previous experiments with bubbles (l3j measured only one component of the velocity and acceleration in a 
von Karman flow. Here, we provide results on the three components of the velocity, acceleration, and autocorrelation 
statistics of micro-bubbles (/3=3) in homogeneous and isotropic turbulence. The micro-bubbles are dispersed in a 
turbulent water tunnel at moderate ReA (160 — 265), and their size is comparable to the Kolmogorov length scale. 
The micro-bubbles can respond to the small-scale fluctuations in the flow and hence the effects of finite-size are not 
important in this study. 

The structure of this paper is as follows: in section [TT] we describe the experimental facility and the smoothing 
algorithm for the particle trajectories. The results are presented in section ITTA followed by a conclusion and summary 
in section HVl 

II. EXPERIMENTS AND DATA ANALYSIS 
A. Experimental Setup 

We conduct experiments in the Twente Water Tunnel, an 8 m long vertical water tunnel designed for studying 
two-phase flows (see Figure [p. By means of an active grid, nearly homogeneous and isotropic turbulence with Re a 
up to 300 is achieved [see Il6l. Il7|. A measurement section with dimensions 2x0.45x0.45 m 3 with three glass walls 
provides optical access for the three-dimensional PTV system. Micro-bubbles with a mean diameter c4 = 340 ± 120 
/j,m are generated by blowing pressurized air through a ceramic porous plate that is located in the upper part of the 
water tunnel. These micro-bubbles are advected downwards by the flow and pass through the measurement section. 
Our three-dimensional PTV system consists of four Photron PCI-1024 high-speed cameras that are synchronized with 
a high-energy (100 W), high-repetition rate (up to 10 kHz) Litron laser (LDY303HE). The four cameras are focused 
at the center of the test section on a 40 x 40 x 40 mm 3 measurement volume that is illuminated by expanding the laser 
beam with volume optics. The arrangement of the cameras and laser is such that the four cameras receive forward 
scattered light from the micro-bubbles. We acquire images at 10,000 frames per second (fps) with a resolution of 
256 x 256 pixels, resulting in a spatial resolution of about 156 \xm/ pixel. 

The ReA is varied by changing the mean flow speed of water in the tunnel. Table U summarizes the flow properties for 
the various cases considered. The flow was characterized by measurements using a cylindrical hot-film probe (Dantec 
Rll) with a sampling rate of 10 kHz placed in the center of the imaged measurement volume. The dissipation rate 
e was obtained from the Kolmogorov scaling for the second-order longitudinal structure function Dll = C2(er) 2 / 3 , 
with C2 = 2.13 [l8j]. For each case of mean flow speed in the water tunnel, the dissipation rate is obtained from the 
value of the plateau region (see Figure [2]), and other parameters follow. 

For the three-dimensional particle tracking, we use the open source code developed at the IfU-ETH group [191 ] . The 
error in the determination of the particle's position is within sub-pixel accuracy of 60 |xm, corresponding to the 
tolerance of epipolar matching in three dimensions. In this paper we focus on the Lagrangian statistics of micro- 
bubbles, but it is also possible to study particle clustering using the three-dimensional data (see [13]). Here, the raw 
particle trajectory is smoothed out with a polynomial fitting method (see section Hi B[) . Velocities and accelerations 
are obtained by differentiating the particle positions in the filtered trajectory. For the Lagrangian statistics shown in 
the results, the number of data points (Ndata) used are larger than 4.5 x 10 6 . 

B. Smoothing method for particle trajectories 

Experimental errors in the determination of the particle positions are unavoidable, and obtaining the particles' 
velocity and acceleration through a time differentiation of their positions would be very sensitive to these experimental 
errors. Hence, a smoothing of the particle trajectory has to be carried out. This smoothing process is a trade-off 
between filtering out the experimental noise and retaining the turbulent features of the particle motion. Therefore, 
the smoothing parameters must be very carefully selected. 
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FIG. 1: The Twente Water Tunnel: an experimental facility for studying two-phase turbulent flows. The picture shows the 
measurement section and on top the active grid, which allows homogeneous and isotropic turbulent flows upto Re^ = 300, 
and the 4-camera particle tracking velocimetry (PTV) system to detect the positions of particles in three-dimensions. For 
illumination we use a high energy, high-repetition rate laser. Micro-bubbles with a diameter «340 urn are generated above the 
active grid using a ceramic porous plate and are advected downwards into the measurement volume. 



TABLE I: Summary of the flow parameters. Vmean'- water mean flow speed, Ke\ = (15i4 ms /e^) 1/ ' 2 : Taylor-Reynolds number, 
u rma : mean velocity fluctuation, 77 = (i/ 3 /e) 1 ^ 4 and r v = (u/e) 1 ' 2 ; are the Kolmogorov's length scale and time scale respectively, 
L: integral length scale of the flow, e: mean energy dissipation rate, St — t p /t v : Stokes number, and Ndata'- number of data 
points used to calculate the Lagrangian statistics. 
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In the Lagrangian Particle Tracking community, two different methods have been used to smoothen particle tra- 
jectories in turbulent flows. One method consists of fitting the trajectory to a polynomial of second or higher order 
IE [HI: whereas the other method uses a Gaussian kernel 0,0, H2|- We tested both smoothing methods, and obtained 
very similar results. In the present paper, we only show results obtained by smoothing the micro-bubble trajectories 
with a third-order polynomial (also referred as the moving cubic spline method). 

The entire signal of the trajectory x(t) is low-pass filtered by fitting a third-order polynomial. Using a fitting 
window with the particle positions from t — Ndt until t + Ndt, where dt is the timestep, the filtered particle position 
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FIG. 2: Flow parameters characterization from hot-wire measurements. Compensated second-order longitudinal structure 
function Dll{t) calculated in order to estimate dissipation. 



at time i is calculated as follows: 

Xij(t) = dfl + C iA t + Ci,2* 2 + Ci^t 3 . 

The Lagrangian velocity and acceleration are obtained by differentiating the particle trajectory: 

(*) = °iA + 2ci,2* + 3c ii3 t 2 , 



(1) 
(2) 



Oi,f(t) = 2cj i2 + 6c ii3 i. (3) 

The parameter N determines the length of the time window (t — Ndt , t + Ndt), and has to be appropriately chosen 
to ensure that the time fitting window is smaller than the typical turbulent time scale. We explore the effect of N 
on the r.m.s (root mean square) of the micro-bubble velocity (Figure [3^) and acceleration (Figure [3}d) to find the 
optimum value for the case of Re,\=195. One can observe in Figure [3b that the r.m.s of the velocity saturates at 
around N = 40, for smaller values of N there is an exponential rise owing to the noise. Since the acceleration is a 
second-order derived quantity, it is more sensitive to the choice of N as observed in figure [2t>. Here, we can clearly 
distinguish two different regions: for small values of N < 30 the r.m.s again increases exponentially due to the noise, 
while at large values of N > 100 the a rms reduces considerably as an effect of the over-smoothing. For the data 
presented in this work we have chosen values of N in the range 45 — 50, which correspond to normalized values of 
iV/fps x r ?; in the range 0.03 — 0.14. It is important to point out that the normalized acceleration PDFs obtained by 
choosing N in this range (45 — 50) are similar for each Re A- The flatness of the acceleration PDF can also provide a 
measure to identify the optimal value of N (as shown in Figure [3b). We describe the details of the flatness calculation 
procedure in section IIII CI Here, in Figure [3J;, we see that our chosen optimal value of N = 50 corresponds to the 
starting point in a region where the flatness values are decreasing with N as a result of over-smoothing. At N < 10, 
the decrease in flatness is an artificial effect arising from the noise. Hence, the optimal value of N is chosen such that 
we do not over-smoothen the micro-bubble trajectories. 



III. RESULTS 



A. PDFs of micro-bubble velocity 



In this section we present results on the PDFs of micro-bubble velocity obtained by smoothing the raw trajectories. 
Figure|4ji shows the PDF of the three components of the normalized micro-bubble velocity at Rca = 195. We observe 
that the velocity distributions of the three components closely follow a Gaussian profile. The flatness values F of these 
PDFs for different Re> (see Table HI)) are close to that of a Gaussian distribution (F = 3). Gaussian- type flatness 
values have also been reported for neutrally buoyant particles in turbulent von Karman flows. Voth et al. |l| measured 
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FIG. 3: r.m.s of the vertical component of the micro-bubble (a) velocity and (b) acceleration, and (c) the acceleration flatness 
at ReA = 195 as a function of N for polynomial smoothing. The arrows in the figures indicate the chosen value (N = 50) for the 
smoothing of the trajectories. 

TABLE II: flatness values of the distribution of micro-bubble velocities. 
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velocity distributions close to Gaussian with flatness values in the range 2.8—3.2, and more recently, Volk et al. [22| 
obtained sub-Gaussian distributions with flatness around 2.4—2.6. In Figure^, we show a plot of the r.m.s values of 
the micro-bubble velocity versus ReA- We observe an increasing trend with ReA, and the three components are nearly 
isotropic. We also compare the r.m.s. values obtained from the hot-wire data and find a reasonable agreement with 
the 3D-PTV velocity measurements. 

B. PDFs of micro-bubble acceleration 

Contrary to the velocity PDFs, the micro-bubble acceleration PDFs normalized with the r.m.s (a/a rms ) exhibit a 
strong non-Gaussian behavior. Figure|4jD shows the PDFs for all the three components of the micro-bubble acceleration 
at ReA = 195. We observe that the acceleration PDFs are highly intermittent with stretched tails that extend beyond 
5 a rms , indicating that the probability of rare high acceleration events is much higher than for a Gaussian distribution. 
At this Rca(= 195), the acceleration is nearly isotropic as the PDFs of the three components show a good collapse 
atleast till 5& rms . We have observed the same trend for the higher ReA, whereas for smaller ReA, the components of 
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FIG. 4: (a) PDFs of the three components of micro-bubble velocity at ReA = 195. The three velocity component distributions 
are nearly Gaussian compared to the solid line that represents a Gaussian distribution, (b) PDFs of the three components of 
the normalized micro-bubble acceleration at ReA = 195. The three components of the acceleration are strongly non-Gaussian, 
i.e. the tails of the distribution show high intermittency. 

the acceleration (a x , a, y ) in the plane perpendicular to the mean flow direction are not yet isotropic. These components 
have tails that are slightly narrower than the vertical component (a z ). The flow in the Twente Water Tunnel is not 
fully isotropic, as has been discussed in Poorte & Biesheuvel (23|. This slight anisotropy is visible in the PDFs, also 
in the a x component in Figure 0Jd. In Figure 03, we show the dependence of the acceleration r.m.s values on ReA- 
Again, there is an increasing trend with ReA and a visible anisotropy in the three components. In the discussion that 
follows, we will only present results of the vertical component z (mean flow direction) of the acceleration PDF. 

In Figure [6l we present the PDFs of the micro-bubble acceleration for all the ReA covered in the present study. In 
order to improve the statistics of the data, we take the absolute value of the acceleration vector. Although the r.m.s 
of acceleration increases with ReA (Figure [SJs), the acceleration PDFs, normalized by the r.m.s, collapse on top of each 
other for all ReA (see Figure [5]). Here, we cannot see a clear dependence on ReA, but the flatness of these acceleration 
PDFs better reveals the dependence, and will be discussed in section UlI CI 

First, we compare our micro-bubble results with experimental data of heavy Q and tracer particles [|| under similar 
flow conditions (grid-generated turbulence with a mean flow). Ayyalasomayajula et al. 0] conducted experiments with 
heavy particles (water droplets) in a wind tunnel at St = 0.15 and ReA = 250, the corresponding results are shown 
in figure [6^,. Subsequently, measurements of tracer particles were also made in the same facility with St = 0.01 and 
ReA = 250 § (also shown in figure^). We observe that the present micro-bubble acceleration PDF shows a higher 
intermittency than heavy and tracer particles at similar ReA and St. 
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FIG. 5: (a) r.m.s values of the three components of the micro-bubble velocity for all ReA, compared with the hot-wire probe 
measurements, (b) r.m.s values of the three components of the micro-bubble acceleration for all ReA- 



In figureinbi we compare the present data with experiments carried out in von Karman flows. It is important to note the 
differences in the flow conditions between grid-generated turbulence and the turbulence generated in between counter- 
rotating disks. It is known that von Karman flows have a large-scale anisotropy. Secondly, the confinement conditions 
are different. These differences could affect the Lagrangian dynamics (Tlj . Tracer particles in von Karman flows at 
ReA=140 to 690 [lj, ||f show a good agreement with our micro-bubble results for \a z /a z , r ms\ ^ 15 (see the Mordant 
et al. 6] fit in figure |6]b) . But beyond this value our micro-bubble acceleration PDFs are slightly less intermittent. 
Furthermore, a comparison with the experiments of Volk et al. [l3j who measured micro-bubble acceleration in von 
Karman flow at Re>, — 850 and St=1.85 is presented in figure [BJd. Unexpectedly, there is a good agreement between 
the two experiments with very different Re A , St, and flow conditions. This could just be a coincidence that both the 
results agree despite these differences. To arrive at a final conclusion on this issue, more systematic experiments need 
to be carried out in a wider ReA- St parameter space under our experimental conditions. 

Figure(SJ: shows the comparison between the present micro-bubble acceleration results with the DNS data for point-like 
bubbles and tracers in homogeneous and isotropic turbulence at a similar ReA=180 and St = 0.1 (data obtained from 
iCFDdatabase http://cfd.cineca.it). The simulations considered one-way coupled point particles within a periodic 
cubic box of size L — 2n and with a spatial resolution of 512 3 (for further details on the simulation see 15]). We 
observe that our experimental findings agree with both numerical bubbles and tracers when \a z /a z ^ rms \ < 15 within 
experimental error. The experimental PDF is closer to the numerical tracers for \a z /a z rma \ > 15. A possible reason 
for the better agreement between experimental micro-bubbles and DNS tracers could be the small St numbers 0(0.01) 
in the present study. Another possible reason is the different flow conditions, in the experiments there is a strong 
mean flow which is absent in the numerics. In addition, several factors such as the lift force, buoyancy forces, and 
particle-particle interactions are ignored in the DNS [24j. We emphasize again that in the present work, the micro- 
bubble size is comparable to Komogorov scale, hence, we do not expect the finite size effects [25j to play an important 
role. 



C. Flatness of the micro-bubble acceleration 



In order to quantify the intermittency of the acceleration PDFs, statistical convergence of the data is necessary. 
The number of data points needed for this convergence is crucial, and previous studies have shown that it should at 
least be re O(10 6 ) |l|, |4j. As shown in table HI our measurements consist of at least 4.5 x 10 6 datapoints. 

The intermittency of the PDFs of the micro-bubble acceleration can be quantified by studying the flatness F: 

flatness = ^ 4 / a A , (4) 

where ^4 is the fourth moment and a the r.m.s of the distribution. The flatness being a fourth order moment is 
strongly determined by the tails of the distribution, and hence convergence of the PDFs is required. Even though 
the number of datapoints used to calculate the PDFs in the present work is larger than O(10 6 ), full convergence has 
not yet been achieved to calculate directly the flatness from the distribution itself (the largest experimental datasets 
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FIG. 6: PDFs of the vertical component of the normalized micro-bubble acceleration, (a) Comparison with experiments under 
similar flow conditions (grid-generated turbulence in a wind-tunnel) at ReA = 250. Our results are shown with open symbols; 
stars are heavy particles 3 and black crosses represent tracer particles Q| . (b) Comparison with von Karman flow results: fit 
for tracers at ReA = 140 — 690 P, [f| is the black line; bubbles at ReA = 850 [lH . fl3l ) are shown with a blue line, (c) Comparison 
with DNS simulations for point particles at ReA = 180 (from iCFDdatabase http://cfd.cineca.it): the red line indicates tracers 
and the black line bubbles. 
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FIG. 7: (a) PDF of the vertical component of the micro-bubble acceleration at ReA = 195. Open squares are the experimental 
data, solid line is the fitted stretched exponential function. The insert shows the plot of the fourth order moment x 4 P(x) for 
experimental data and fit. (b) The flatness value of the fitted PDFs of micro-bubble acceleration as a function of Ke\. (c) The 
flatness values versus the Reynolds number. Comparison with Voth et al. [IJ and Ishihara et al. [§[ reveals that the present 
micro-bubble trend agrees well with the data in the literature, at least till Re\ = 225. 



consist of w O(10 8 ) datapoints @, whereas for numerics this value can go up to ~ O(10 9 ) Q. Consequently, we fit 
the experimental PDF to a stretched exponential distribution [6, 26] defined as: 

f(x) = Cexp( |, (5) 

V« 2 ( 1 + lfl 7 )/ 

In equation [5l x — a/a rms is the normalized micro-bubble acceleration, the fitting parameters are a, f3 and 7 and C 
is a normalization constant. For this fitting procedure and in order to improve the convergence at the tails, we have 
taken the absolute value of the acceleration vector. 

Figure [JJi shows the result of the fitting for the micro-bubble acceleration PDF at Re A = 195. The stretched 
exponential fits the experimental PDF quite well because the three fitting parameters enable a fine adjustment. In 
the inset of figure [7^,, we plot the fourth order moment (a/a rms ) 4 PDF(a/a rms ) for the experimental acceleration 
measurement along with the fitted curve. This type of curve allows for a good convergence test (27j . At the tails of 
the distribution, convergence is nearly achieved, and the fitted curve nicely sits on top of the experimental data. We 
have observed a similar behavior for the other measurements at different Rc A . 

Next, we calculate the flatness of the fitted acceleration PDFs as a function of Re A , as shown in Figure [JJo. The 
flatness is determined directly from the fitted stretched exponential functional for all Re a- The errorbars are obtained 
by finding the difference between the flatness values for half and the entire acceleration datapoints. Figure [7J: shows 
that the flatness values of the micro-bubble acceleration PDFs increase in the Re A range 160-225 consistent with 
the experimental results of Voth et al. [l[ for tracer particles and the numerical results of Ishihara et al. [|| for fluid 
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particles. For the highest ReA, we have less statistics compared to the other cases as the mean flow speed is the 
fastest. This might be the reason for the decrease (underestimation) in the flatness value of our data point at ReA = 
265. Clearly, from the collection of all data (Figure [Th) one would not expect such a decrease. 




FIG. 8: Autocorrelation function of the three components of the micro-bubble acceleration at ReA = 195. The acceleration 
autocorrelation of the micro-bubbles is nearly isotropic. The time lag is normalized with the Kolmogorov time scale t v . 



D. Autocorrelation functions 

We now present results on the Lagrangian autocorrelation function of the micro- bubble acceleration. In figure [5] we 
compare the autocorrelation for the three components of the acceleration at ReA=195, using a time lag normalized 
with t,,. We find that the three acceleration components correlate in a similar manner. This nearly isotropic behavior 
was also found for the other measurements at different ReA- 

Figure^ shows the autocorrelation of a z for different i?eA- It is clear that the microbubble's acceleration correlates 
for longer times as ReA increases. The acceleration autocorrelation function drops to zero rapidly and the zero- 
crossing point has small values: < O.It^. Voth et al. [l[ and Mordant et al. 0] reported values of around 2.2^ in their 
experiments with tracers in von Karman flows at high turbulence intensities (ReA > 690). The value of 2.2^ was first 
found from DNS by Yeung et al. (28J. Volk et al. [l2j performed both micro-bubble and tracer experiments in a von 
Karman apparatus, and found that the decorrelation of the microbubbles is smaller than that of tracers at a given 
ReA- We do not yet know the exact reason for the large disparity between 2.2^ for the fluid particles compared to 
O.Itv, for the present micro-bubbles. One possible reason is that our flow conditions are different as we have a strong 
mean flow. 

We study the time at which the autocorrelation function drops to zero for different ReA by defining the decorrelation 
time as: 

T D = [ ° C a (T)dT, With C a {To) = 0, 

Jo 

where C a is the acceleration autocorrelation function. Tjj represents the characteristic time for the evolution of the 
micro-bubble response to changes in the flow conditions. FigureEb shows the dependence of Td /t v on ReA for the three 
components of the micro-bubble acceleration. We observe that Tjj /t v increases with ReA , and that the autocorrelation 
functions are nearly isotropic as evidenced by the very similar Tjj / values for the different components. In the inset 
of figure [9b, the decorrelation time Td/t v as obtained by Volk et al. [l2j at ReA=850 agrees well with our increasing 
trend of the decorrelation time with ReA ■ We fit our experimental data of the decorrelation time of a z to a linear 
relation Td/t, ; = 0.00038Rca — 0.051 that is shown in the inset of figure [Hb as a solid line. Evaluating these relations 
at ReA=850 gives Td/t v = 0.27, which agrees well with their experimental value of Td/t v = 0.258 [l2|. More 
experiments are needed to fill the gap of ReA ■ Very recently, Volk et al. [22j found an increase of Tu /r„ with ReA for 
a fixed particle size, just as we find in our micro-bubble experiments. 
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FIG. 9: (a) Autocorrelation function of the vertical component of the micro-bubble acceleration for the different ReA measured. 
The correlation of the micro-bubble acceleration persists longer with increasing Reynolds number, (b) The decorrelation time 
Td/t v of the autocorrelation function for the three components of the micro-bubble acceleration as a function of ReA- The 
decorrelation time increases with the turbulent intensity. In the inset, we show also the result of Volk et al. il2j at a very high 
ReA=850 (♦), their experimental point agrees with the trend of increasing decorrelation time with turbulent intensity. The 
linear fit obtained with our experimental data extrapolates a value of Td/t v = 0.27 at ReA=850, which is slightly higher than 
their experimental value of Td /t v = 0.25. 



IV. CONCLUSION 



We have presented experimental results on the Lagrangian statistics of micro-bubble velocity and acceleration in 
homogeneous isotropic turbulence. Three-dimensional PTV was employed to obtain the micro-bubble trajectories. 
The PDFs of micro-bubble velocity closely follow a Gaussian distribution with flatness F ~ 3, independent of ReA- But 
the acceleration PDFs are highly non-Gaussian with intermittent tails. Although the acceleration PDFs themselves 
do not show a clear dependence on ReA, the flatness values reveal a clear trend. We fit the experimental acceleration 
PDFs to a stretched exponential function and estimate the flatness based on the fitting. The flatness values were 
found to be in the range of 23-30 and show an increasing trend with ReA- This trend is consistent with previous 
experimental [l[ and numerical Q results. 

A comparison of our results with experiments in von Karman flows [l], 0, Il2l ll3T | suggest that the present micro- 
bubble acceleration PDF is similar to tracers and bubbles (in von Karman flows) for very different Re\. However, 
there are significant differences in the flow conditions between the two experimental systems. Therefore, it is more 
relevant to compare our results with previous investigations in similar flow conditions, i.e. grid-generated turbulence. 
We find that the acceleration PDFs of our micro-bubbles are more intermittent as compared to heavy and tracer 
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particles in wind tunnel experiments at similar St and Re a [3|, |j| . 

Compared to DNS simulations in the point particle limit, our micro-bubble acceleration PDFs show a reasonable 
agreement with both numerical tracers and bubbles, but in the tails our data has a better match with numerical 
tracers. One possible reason is the differences in flow conditions between the experiments and numerics. Another 
possibility is that the St in our experiments are small (0.02 — 0.09). It will be interesting to study the acceleration 
statistics of finite-sized bubbles at large St in turbulent flows. 

We also calculate the autocorrelation function of the micro-bubble acceleration, and observed that the decorrelation 
time increases with Re>. This finding is consistent with other experimental investigations [12| . |22| at very high Reynolds 
number. More experimental data is needed to fill the gap of Re> in order to further study the scaling behavior. 
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